Self-consistent Eliashberg theory, T c , and the gap function in electron-doped cuprates 
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We consider normal state properties, the pairing instability temperature, and the structure of 
the pairing gap in electron-doped cuprates. We assume that the pairing is mediated by collective 
spin excitations, with antiferromagnetism emerging with the appearance of hot spots. We use a 
low-energy spin-fermion model and Eliashberg theory up to two-loop order. We justify ignoring 
vertex corrections by extending the model to N >> f fermionic flavors, with 1/N playing the role of 
a small Eliashberg parameter. We argue, however, that it is still necessary to solve coupled integral 
equations for the frequency dependent fermionic and bosonic self-energies, both in the normal and 
superconducting state. Using the solution of the coupled equations, we find an onset of d— wave 
pairing at T c ~ 30 K, roughly three times larger than the one obtained previously [P. Krotkov 
and A. Chubukov, Phys. Rev. B 74, 014509 (2006)], where it was assumed that the equations for 
fermionic and bosonic self-energies decouple in the normal state. To obtain the momentum and 
frequency dependent d-wave superconducting gap, A(kp, w„), we derive and solve the non-linear 
gap equation together with the modified equation for the bosonic self energy which below T c also 
depends on A(kF,L} n )- We find that A(fci?,w n ) is a non-monotonic function of momentum along 
the Fermi surface, with its node along the zone diagonal and its maximum some distance away from 
it. We obtain 2A max (T — > 0)/T c ~ 4. We argue that the value of T c , the non-monotonicity of the 
gap, and 2A max /T c ratio are all in good agreement with the experimental data on electron-doped 
cuprates. 

PACS numbers: 



I. INTRODUCTION 

The first high temperature superconductors, discov- 
ered by Bednorz and Miiller in 1986 1 -, were hole-doped 
cuprates, with Sr doped into La2Cu04. Electron-doped 
cuprates were subsequently discovered in 1989, by dop- 
ing Ce into Nd2Cu04— . At present, there exist a variety 
of electron-doped cuprates of the form RE 2 _ a; M :l; Cu04 
where RE (rare earth lanthanides) is Nd or Pr and M is 
Ce or Th. For a comprehensive review on electron-doped 
cuprates see RefA 

Raman spectroscopy 4 ^ and ARPES 6 '"- 9 - provide strong 
evidence that the symmetry of the superconducting 
gap in both electron- and hole-doped cuprates is 
d x 2_ y 2. This gap symmetry is also consistent with 
tunnelling 10 , penetration depth 1 1-14 , and Andreev re- 
flections measurements^, although some of these exper- 
iments were originally interpreted differently On the 
other hand, maximum T c in electron-doped cuprates is 
in the range of 10 — 30 K, nearly an order of magnitude 
smaller than in hole-doped cuprates despite the "Hub- 
bard U" being the same in both sets of materials (optical 
conductivity studies of parent hole- and electron-doped 
materials reported the same value, 1.7 eV, of the optical 
gap 16 ' 17 ). Additionally, the d-wave gap in at least some 
electron-doped materials shows a non-monotonic behav- 
ior along the Fermi surface (FS), with the maximum 
somewhere in between the zone diagonal (where the gap 
vanishes) and the region near (0,tt) (Fig. [TJ. This non- 
monotonicity of the d— wave gap was predicted based on 
the analysis of Raman dataA, and was subsequently mea- 
sured directly by ARPES in Pr .89LaCeo.iiCu0 4 (Ref. 9 ). 

One way to obtain d x 2_ y 2 pairing is to assume that the 



pairing is of electronic origin and is mediated by collective 
bosonic excitations at large momentum transfer, such 
that the interaction predominantly couples fermions with 
gaps of different signs (this effectively converts repul- 
sive interaction into an attractive one) . Both charge and 
spin fluctuations can give rise to such pairing. The idea 
of spin-fluctuation mediated pairing in cuprates has at- 
tracted substantial attention because antiferromagnetism 
is part of the phase diagram for both hole- and electron- 
doped cuprates. We follow earlier works on both hole and 
electron-doped cuprates 1 -^— and other materials 2 ^—, 
and consider spin-mediated pairing. 

In this paper we consider two issues. First, the origin 
of the smallness of T c in electron-doped cuprates com- 
pared to hole-doped cuprates; second, the momentum 
and frequency dependence of the superconducting gap. 

The relative smallness of T c in electron-doped cuprates 
is generally attributed to the fact that in these systems 
d— wave pairing is less effective because hot spots (points 
at the FS separated by the antiferromagnetic momen- 
tum Q = (jt,tt)) are located much closer to the zone 
diagonals (see Fig. [T|). This argument, however, is in- 
complete because magnetically-mediated T c remains fi- 
nite and scales linearly with the overall spin-fermion cou- 
pling, even when hot spots merge at the zone diagonals 28 . 
The complete picture is that in electron-doped cuprates, 
T c acquires the dependence of the FS curvature because 
the velocities of hot fermions become nearly antiparallel 
to each other when hot spots approach zone diagonals 
(Fig. [1} . Krotkov and one of us 2 ^ (hereafter referred to 
as KC) demonstrated that it is this dependence that gives 
rise to the additional numerical smallness of T c . Using 
the same Eliashberg-type computational procedure as in 
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FIG. 1: Brillouin zone for electron-doped cuprates at the dop- 
ing when hot spots merge on the zone diagonals. We argue 
that this doping is close to the one at which antiferromagnetic 
order emerges. The arcs represent the FSs, the dashed square 
represents the antiferromagnetic Brillouin zone. k x and k y 
are the directions of fermionic momenta transverse and along 
the FS, respectively. 



hole-doped cuprates, KC obtained T c ~ 10 K, using for 
the coupling the value extracted from the optical data at 
half-filling. 

Compared to experiment, the value of 10 K is a bit 
small, particularly given that "hot spot constrained" cal- 
culations tend to somewhat overestimate T^. 

In light of these points, we re-consider here the Eliash- 
berg computational procedure. We argue that the ear- 
lier study neglected corrections of order one and should 
be therefore modified. To understand where the mod- 
ification is required, we note that for superconductiv- 
ity mediated by collective boson exchange Eliashberg- 
type calculations of T c differ qualitatively from those 
of phonon-mediated superconductivity. For phonon su- 
perconductivity, the full phonon propagator is an input 
with the bosonic self-energy neglected because it is small 
by the same parameter as vertex corrections. For col- 
lective mode mediated pairing, bosonic dynamics origi- 
nates from low-energy fermions and has to be computed 
together with the fermionic self-energy. As a conse- 
quence, one generally has to solve a coupled set of integral 
equations for the fermionic self-energy, E(fc,w), and the 
bosonic self-energy (the bosonic polarization operator), 

For hole-doped cuprates, the solution of this coupled 
set near hot spots is further simplified because bosonic 
self-energy, H(Q,fl), has the form of Landau damping 
and does not depend on the fermionic self-energy. As a 
result, the set of coupled equations for E(w) = E(fch s ,w) 
and n(<5, £1) reduces to just one equation for E(w), while 
Tl(Q,Q) remains the same as for free fermions. 

The reasoning for the independence of IT(Q, from 
E(w) is the following: n(Q, O) is aconvolution of a par- 
ticle and a hole, located near hot spots separated by 
Q. When these two fermions have Fermi velocities di- 
rected at an arbitary angle w.r.t. each other, each 
fermionic line can be integrated over its dispersion 
ffc+Q, and Tl(Q, f2) emerges as the convolution of the two 



local propagators Q\{lo) = f dei~G(k, cj) and Qi(uj + fi) = 
/ dek+QG{k + Q,uj + ft). The local Q\{uj) — — i7rsgn(w), 
and does not depend on fermionic E(cj); hence one does 
not need to dress up fermionic propagators. An alternate 
way to state this is to observe that Landau damping is an 
anomaly, i.e., it can be re-expressed as the contribution 
from high-energies, where fermions are free particles. 

Applying Eliashberg approach to electron-doped 
cuprates, KC assumed that the equations for E and IT de- 
couple in this case too, so that they used the free-fermion 
H(Q,ui) to compute E(w) and T c . However, H(Q,ui) for 
electron-doped cuprates scales as y]f2[ rather than as 
|Q|. This is because the Fermi velocities of hot fermions 
are anti-parallel, and the Jacobian of the transformation 
from J d 2 fc to J de^defc+Q diverges, so that n(Q,f2) no 
longer reduces to the convolution of local propagators 
(alternatively, H(Q,il) is not an anomaly). As a conse- 
quence, H(Q,Sl) acquires the dependence on E(w), and 
one has to solve the coupled set of equations for E and 
IT. 

We found that to a good numerical accuracy, the so- 
lution of the coupled set modifies previous results of KC 
by a constant factor a ~ 0.6. That is, H(Q, ft) by a and 
E gets multiplied by 1/ \fa. Recalculating T c using these 
rescaled quantities, we find that T c increases by a factor 
1/a 2 ~ 3. This brings Eliashberg T c to 30 K, which is 
quite consistent with the experimental data. 

This T c ~ 30 K is also in better agreement with T c 
obtained within the FLEX approximation^. In FLEX, 
full Green's functions are used, but vertex corrections are 
neglected. The fact that for electron-doped cuprates we 
also need to use the full Green's functions makes Eliash- 
berg theory more of a match with FLEX than for hole- 
doped cuprates. Nevertheless, Eliashberg and FLEX are 
still not equivalent because in FLEX the full momentum 
dependent E(fe, ui) is used (and vertex corrections are ne- 
glected, though without any parametrical justification), 
while in our Eliashberg theory, only E(fci?,w) is relevant, 
while dE(fc, w)/defe and vertex corrections are neglected, 
by analogy with the electron-phonon problem. There is, 
however, no "natural" Eliashberg parameter analogous to 
ujd/ Ep, and to justify the approximation we need to ex- 
tend the theory to N fermionic flavors and consider large 
N limit; then 1/N becomes the Eliashberg parameter. 

Very recent studies^ have demonstrated that the large 
N approximation in fact does not work beyond two-loop 
order. There are a series of multi-loop diagrams, begin- 
ning at three-loop order, which do not contain 1/N, al- 
though apparently small numerically. We restrict to two- 
loop order, and do not address this issue in the present 
MS. 

The second issue we consider is the form of the su- 
perconducting gap as a function of momentum along the 
Fermi surface. KC obtained a non-monotonic gap func- 
tion by solving for the <i-wave eigenfunction of the gap 
equation at T c . Here we solve the non-linear gap equa- 
tion for T < T c (together with the coupled equation for 
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the bosonic H(Q,Cl)), and find that the gap remains non- 
monotonic for T < T c ; there is a slight change in shape 
with T, but the non-monotonic character remains intact. 
At T — > 0, we find that the maximum value of the gap 
along the FS, A max , is about two times larger than T c , 
i.e., within Eliashberg theory, 2A max /T c ~ 4. This ratio 
is quite consistent with these experimental data: from 
optics we have-21 2A max /T c ~ 5 for Pri.gsCeo.isCuO^ 
and from Raman measurements^ on P^-zCe^CuC^ and 
Nd 2 _a;Ce x Cu04 we have 2A max /T c ~ 3.5. 

The paper is organized as follows. In Sec. II we 
briefly review the spin-fcrmion model for electron-doped 
cuprates and its extention to N » 1. In Sec. Ill we 
carry out normal state analysis. We show that vertex 
corrections are small by 1/N, but that self-energy correc- 
tions to n(Q, f2) are important. We then solve coupled 
equations for E(cd) and H(Q, f2), and compare our results 
with those of KC, who used bare fermionic propagators 
to compute H(Q, f2). We find that the fermionic self- 
energy scales as {Cd ) 1 / i \uj\ i / i , and we use this wq scale 
as a measure for T c . In Sec. IV we study the super- 
conducting properties. We first solve for T c , and then 
derive and solve the non-linear gap equation. The latter 
yields the non-monotonic d— wave gap A.(kp,uj n ), which 
depends on momentum along the FS and also on Mat- 
subara frequency. Using this solution, we obtain the ratio 
2A max /X' c . In Sec. V we summarize and discuss our re- 
sults. A vertex correction calculation is relegated to the 
Appendix. 



II. THE MODEL 

Spin-fermion model has been discussed before, so we 
will be brief. The idea is that the low-energy physics of 
a system of itinerant fermions near antiferromagnetic in- 
stability, is adequately described by the interaction with 
low-energy, collective, bosonic excitations in the spin 
channel. The static part of a collective mode propaga- 
tor comes from fermions at high energies and is the in- 
put for the low-energy theory, while the dynamical part 
comes from low-energy fermions and has to be calculated 
within the low-energy theory. The Hamiltonian of the 
spin-fcrmion model is 

H = X! efcC U C fc,a +J2Xst\<l)Sq ■ S -q 
k q 

+ 9 Cfc+g.a^a./SCfc^ ' S- q 



(i) 



where e& is the electronic dispersion, Ck.a is the fermionic 
operator for an electron with momentum k and spin a, 
a are the Pauli spin matrices, and S q is the bosonic vari- 
able describing collective spin degrees of freedom. Fur- 
ther, g is the coupling constant and Xst(l) is the static 
spin susceptibility (the static propagator of collective 
spin degrees of freedom). The triple spin-fermion ver- 
tex is shown diagrammatically in Fig. [5J We will not 




FIG. 2: The vertex for the spin-fermion interaction. The solid 
lines are fermions, the wavy line is a spin fluctuation. 



show the spin indices explicitly in subsequent Feynman 
diagrams; it should be understood that they are there, 
and proper accounting of them has been done to calcu- 
late the various prefactors while translating the diagrams 
into equations. 

For the static spin susceptibility, we follow earlier 
works and use the Ornstein-Zcrnike form 



Xst(q) 



Xo 



r 2 + (q - Q) 2 



(2) 



where £ is the magnetic correlation length. The constant 
X should not be taken as another input because g and 
Xo only appear in the combination g = g 2 Xo- Through- 
out the paper we refer to g as the effective spin-Fermion 
coupling constant. 

Numerical studies of the onset of antiferromagnetic 
instability all show that antiferromagnetism emerges 
around the doping when hot spots merge along zone di- 
agonals and the Fermi surface touches the antiferromag- 
netic Brillouin zone, as shown in Fig. [T] The analytical 
argument is that the derivative of the static susceptibility 
of free fermions XstiQ) with respect to doping diverges 
at this point, so that in its near vicinity UXstiQ) reaches 
one and within RPA the full Xst(Q) oc l/(l-Ux°t(Q)) di- 
verges. We assume for definiteness that the antiferromag- 
netic quantum critical point (QCP), at which =0, is 
right at this doping. Our results will change only a little 
if £ diverges not at this point, but close to it. 

As mentioned in the Introduction, the velocities of the 
fermions along Brillouin zone diagonals are anti-parallel 
to each other, with the tangential component along the 
FS vanishing (see Fig. [1]). The dispersion of fermions 
near the hot spots is then not the usual = vp x k x + 
VFyky, but 



e k = v F k x + P 2 k 2 



e k+Q 



-v F k x 



P 2 k 2 



(3) 



(4) 



where ft 2 parameterizes the curvature of the Fermi sur- 
face; this will play an important role in our analysis. 

It is convenient to measure the combined effect of the 
curvature and the interaction in terms of the dimension- 
less parameter 



(5) 



For parameters relevant to electron-doped cuprates, r ~ 
10 _1 . We also introduce for further use, the momentum 
scale 
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TTVf 



and the frequency scale 

( 3P x 1 



w 



qo/3 

2 



4-7T 



(6) 



(7) 



The computational procedure is straightforward. We 
need to obtain the dynamic spin susceptibility, x(<T, fin), 
in which the frequency dependence comes through the 
total polarization operator n to t(9, fin), which contains 
contributions from both umklapp and non-umklapp pro- 
cesses. We have 



(8) 



(g-Q) 2 + Xon t ot(«,fin) 
where 

n n = 2nirT (9) 
is a Bosonic Matsubara frequency, and 

n tot (<f, fi n ) = 2(n(£ n n ) + u(-q, a n )) (10) 

We also need to obtain the full normal and anomalous 
Green's functions Q{p) and F(jp), respectively. These de- 
pend on the normal and anomalous self-energies, E and 
E02, as 



Q = 
T = 



X(- P ) 



X{-p)X{p) + {zZ m {p)f 
E 02 (p) 

x(-p)x( p ) + (j:o2(p)) 2 



(ii) 



where p = (p,ui n ), ui n = TtT(2n + 1) is a fermionic Mat- 
subara frequency, and 

x( P ) = go 1 (p) + ^(p), g 1 ( P ) = ^n-e p . (12) 

In Eliashberg approximation, which we justify below, 
the normal and anomalous fermionic self-energies are 
given by one-loop diagrams (Fig. [3]) that generally in- 
volve the full fermionic propagators and the full dynam- 
ical spin susceptibility; similarly, the bosonic self-energy, 
H(q, fi n ), is the sum of bubbles made of the full fermionic 
propagators (Fig.[4j. We have, 



£(p) = -3g 2 T^ J -^^g(k,uj m )x(k - p, 



(13) 



v f d 2 k / 

U(q,Q n ) =2g 2 T22 / / 2 \ 2 \ Q{k,u m )Q{k + q,u m +Q n 

(14) 

These equations can be further simplified by integrat- 



P k P 

x 



-p^>p 



-P -k 



k P 



FIG. 3: Diagrams for the normal and anomalous fermionic 
self-energies (E and E02, respectively) in the Eliashberg-type 
theory (no vertex corrections). The double lines are full 
Green's functions of intermediate fermions, with self-energies 
included. 
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FIG. 4: Diagrams for the polarization operator. The diagram 
on the right is present only in the superconducting state. 

ing the r.h.s. over e&. The dependence on e*, is in 
both fermionic and bosonic propagators, but the one in 
the bosonic propagator can be neglected because over- 
damped bosons are slow modes compared to fermions; 
keeping the dependence on in x, only gives rise to 
a small correction in the Eliashberg parameter. Within 
the Eliashberg approximation, we also neglect the de- 
pendence of the fermionic self-energies on . With these 
steps the integration over e/~ is straightforward. For the 
normal state self-energy at a finite temperature we have 



= i3.9 2 r 

E 

\n m \<\ui„\ 



sgn(w„)x 



dPyX(ky,Py,tt m ) 



(15) 



where k y is the momentum component along the 

FS (v F k x = -P 2 kl) and x(ky,p y ,^ m ) = X{k ~ 
p, f2 m )| £fc=e = o is the spin susceptibility at momenta con- 
necting two points on the FS. Similarly, the linearized 
equation for the anomalous self-energy, valid at T = T c , 
is 

4ttv f ^ 



E 2(fey, U>r, 



dp 



X(ky,Py,U m -W„) 



\ui m + lE(Py,LJ m )\ 



E 2(Py, 0J m ) 



(16) 
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We assume d— wave pairing for which T,Q2{k y ,0J v ) = 
— So2(fcj, + Q,u n ). Shifting p y by Q in Eqs. (fT5|) and 
(|16p . we eliminate the overall minus sign on the r.h.s. 
of (fT6]) . The polarization operator becomes, after the 

shift, n to t(fc — P ~ Q, fl) and can be approximated by 
n to t(Q, ft) = 4n(Q, fi). The susceptibility spin suscepti- 
bility, x(ky,p y ,£l), then takes the form 



(k y - Py y + {p/v F y(ki + pif + 4 Xo n(Q, n m ) (17) 

Xo 



f (fc - p) 2 + r 2 (fc 2 + p 2 ) 2 + %n(Q, ^ 



As mentioned in the Introduction, for fermion-fermion 
interaction mediated by a collective mode, there is 
no natural small Eliashberg-type parameter, analogous 
to lod/Ef, that would make vertex corrections and 
dS(/c)/defe much smaller than one. To rigorously justify 
Eliashberg approximation, we need to extend the theory 
to N > 1 fermionic flavors and take the N >> 1 limit. 
Applying this extension to the low-energy model of Eq. 
([1]), we find that the polarization operator acquires an 
overall factor of N. This factor also appears in the nor- 
mal state fermionic self-energy, and can be absorbed into 
the renormalization of the frequency scale loq (given by 
Eq. J?])), by redefining, 



w 



up 
N 2 



(18) 



Throughout the paper we will use this redefined ujq, keep- 
ing it fixed by rescaling the curvature /? by N. 



A. Vertex corrections 

We first show that diagrams with vertex corrections 
are small by 1/N. As an example, consider the vertex 
correction diagram for the polarization operator (Fig.[5]). 
We have 




FIG. 5: The polarization operator with the vertex correction, 
but without self-energy corrections to intermediate fermions. 



n v (Q,n) = -2Ng 4 



dwd 2 fcdw'd 2 fc' 



, (2tt) 3 (2tt) 3 

G(k + Q, uj + n)G(k', u')G(k' + Q, J + n) x (k -k',u)-u') 

(21) 

where the subscript v stands for II with vertex correc- 
tion. The calculation is presented in Appendix A, and 
the result is 



n v (Q,n) 



9 



'-log Uv^W^ 
Airfivp V 8 8 V V ft 



(22) 



Comparing this with the result for Hq(Q, f2) (Eqn (fl9|) ). 
we immediately see that 



III. NORMAL STATE ANALYSIS 

We begin by citing the results from Ref. [28| for the 
normal state polarization operator and the normal state 
fermionic self-energy at the hot spot, obtained without 
self-consistency. The polarization operator is given by 



n (Q,ft„) 



Ng 2 
2n/3v F 



T 



(19) 



where F n (x — > 0) = 1, Fn(x » 1) ~ y/x, and the 
subscript is used to emphasize that IIo is obtained using 
bare fermionic propagators. The fermionic self-energy at 
the hot spot is given by 



E(w„) = iM 1/4 K| 3 / 4 sgnK)F s ^ 



(20) 



where << 1) oc x a l and Fj:(x >> 1) « 1. For the 

pairing problem we need the region \uj n \ > wo, for which 
Fs can be approximated reasonably well by one. 

We now discuss how these results are modified when 
we go beyond bare vertices and free fermion propagators. 
For definiteness, for the rest of this section we restrict our 
analysis to T — >• 0. 



IL, _ log (4V2N 2 ^) 

n " 



4y/2N 



(23) 



i.e., the vertex correction diagram for the polarization 
operator contains additional log N/N and can be safely 
neglected at large N. One can verify that the same is 
true with the vertex correction diagram for the fermionic 
self-energy. 



B. Self-energy corrections to the polarization 
operator II 



k + Q, lj + v 




FIG. 6: The polarization operator without vertex corrections, 
but with full proparators of intermediate fermions. 
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We next show that corrections from inserting self- 
energy into the fcrmionic lines in the particle-hole bubble 
are not small in 1 /N and should be included. For this, we 
evaluate the polarization bubble II SC (Q, f2) at T — > with 
the renormalized Green's functions (Fig. [6]) and compare 
the result with U (Q, We discuss full self-consistency 
later. The Green's function used is, 



G(k,u) n ) = - 



1 



i (w„ + s(w„)sgn(w„)) - e fc 



where 



a(u n ) = (w»,) 1/4 M 3/4 



(25) 



(24) For the diagram for II SC we have 



U sc {Q,n) = 2Ng 2 



dwd 2 fc 



(27r) 3 (efc - i(u + s(w)sgn(a;)))(efc + Q — i(w + fi + s(u + £l)sgn(u; + fi))) 



(26) 



r 



Expanding the dispersion around the hot spots, carrying 
out the integration over momenta, rescaling the resulting 
expression, and subtracting the zero frequency contribu- 
tion (the static part is already absorbed in Xst(Q))> we 
obtain 



n se (Q,r») 



Uo(Q,n)I 



(w ) ' 



(27) 



where 





f r°° 






oo 





dy 



o v/2y + x + j/ 3 / 4 + (x + y) 3 / 4 



(28) 



/o N /2 2 / + 2y 3 / 4 y' 
We plot I(x) in Fig. [71 This function monotonically in- 




FIG. 7: Plot of the function I{x) given by Eqn. vs. x = 
f2/cjo for x > 1. Observe that I{x) changes little between 
x = 1 and x ~ 10. At very large x, I(x) approaches one. 



creases with x and approaches one at large x. This lim- 
iting behavior can be easily extracted from the integral 
for I(x). 

We see that n so (<5, f2) is of the same order as IIo(Q, O), 
i.e., there is no factor \/N between them. This is not 



surprising as for typical Q ~ ujq (i.e., for typical x ~ 0(1) 
in I(x)), to and are of the same order. This is in 

contrast to the vertex correction diagram in which the 
insertion of an additional bosonic propagator adds a 1 /N 
factor. This factor could potentially be compensated by 
additional powers of momentum and frequency in the 
denominator; however it turns out that the integration 
over momentum and frequency in the diagram for II V 
only adds an additional log N. 

We therefore conclude that it is possible to ignore 
vertex corrections; but two-loop self-energy corrections 
to the polarization operator cannot be ignored. There- 
fore, we do need to calculate S(w) and H(Q,tt) self- 
consistently. In this respect, the calculation of E and II 
in electron-doped cuprates is qualitatively different from 
that in hole-doped cuprates, where H(Q,fl) cx |0| does 
not depend on E, i.e., the equations for E and II decou- 
ple. 

We discuss our self-consistent solution in the next sub- 
section. 



C. Self-consistent analysis of E(cj) and I1(Q, SI) 

For the self-consistent calculation of the bosonic self- 
energy, II, and the normal state fermionic self-energy, 
E(w), we use equations (fT5)l . flTTl) . and the diagram in 
Fig. |H1 In general, this would require one to solve the 
resulting coupled integral equations. In our case how- 
ever, there is a way to simplify this calculation because 
the function I(x) in (|28)l is quite flat at x > 1, varying 
only by 8% between x = 1 and x = 5 (roughly from 0.6 
to 0.65). To a reasonable accuracy, we can then approx- 
imate I(x) by a constant a ~ 0.6. Once we make this 
approximation, the polarization operator, II, can be cast 
into the form 



n se (Q,r!) 



4xo 



(29) 
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where 



rv2 



(30) 



In other words, the dressed II se differs from Ilo in Eq. 
(|19[) by the replacement ujq — > Qq, without changing 
<7o. Because ujq — q 2 /3 2 /4, this renormalization implies 
that the curvature (3 gets effectively renormalized into 
/3 = (3 1 a. Substituting this polarization operator into 
equation (|15[) for the self-energy, we obtain, for uj > u>o, 



£(w) 



\— 2 sgn(o;) 



M 



dft 



dk 



i— sgn(w) 



i( w o) 



1/4 



I" 



3/4 a 



dfi 



3 (fc 2 

d/c 



4xo 



n x (Q,n) 



sgn(u;) 



(31) 



Just as with II SC , the normal state fermionic self-energy 
preserves the form of Eq. (|20[) . with wo replaced by wo- 
Thus, by approximating I(x) as a constant (= a), the 
full solution of the set of self-consistent equations for £ 
and II(= Il sc ) reduces to the replacement ujq — > Wo, with 
no change in go- 



IV. SUPERCONDUCTING PROPERTIES 

A. The value of T c 

The linearized gap equation (equation (([16])) without 
self-consistent renormalization of £ and II (i.e., with ujq 
instead of u>o) has been solved by KC (Refi^£). These 
authors obtained 



T. 



KC 



rg_ 

-in 



F(r) 



(32) 



where, we remind, r — gj3 2 / '(nVp) is the dimensionless 
parameter proportional to the FS curvature. The func- 
tion F(r) is a non-monotonic function of r, with a flat 
maximum near r ~ 10" 1 . Using g ~ 1.6 eV (extracted 
from the charge transfer gap) and t — t' dispersion with 
parameters taken from ARPES, KC obtained ujq ~ 10 
meV and r - 0.08, with F(r) - 0.1 (Ref^). This yields 
T C KC - 10 K. 

As we just discussed, the self-consistency of E and II 
changes the scale ujq into ujq — loq/o 2 , with a ~ 0.6. The 
static part of the susceptibility in (|17p . which provides 
the r— contribution in F(r), is not affected by this renor- 
malization, i.e., F(r) does not change. We then obtain, 



T c (r) = cd F(r) 



T KC 



3T 



KC' 



so that T c changes by a factor of about 3. For the 
same parameters as used by KC, we obtain T c ~ 30 



Tc/g 




0.03 0.04 0.05 0.06 0.07 0.08 0.09 



FIG. 8: T c in units of g versus r within the constant a ap- 
proximation. The analytical formula is T c /g — rF(r)/(4iva 2 ), 
where F(r) is a complicated function of r (see Ref— ). The 
solid line is for M = 20, the dashed for M = 100, where 
M is the number of Matsubara points. The two curves are 
in good agreement, thus providing a check on our solution. 
We kept the number of points for the Gauss-Legendre (GL) 
quadrature fixed at Mgl = 10. 




FIG. 9: Y, m {k y 



versus k y /qo, obtained as a normalized 



eigenvector of the linearized gap equation. We set r — 0.08 
(T c ~ O.SOuj). The magnitude of this £02 does not have any 
significance as the actual £02 is infinitesimally small at T c . 
The n in the legend correspond to different Matsubara fre- 
quencies, CJ n . 



K. For completeness we also computed T c (r) and the 
eigenfunction A(k yi uj n ) numerically, within our constant 
a approximation. We used the Gauss-Legendre quadra- 
ture method, one of the many Gaussian type quadrature 
(33) approximations—. We present the result for T c (r) in Fig. 
[5] We find that the functional form of T c (r) is quite simi- 
lar to the one obtained by KC, and the values are about 3 
times larger than theirs, as expected. We also emphasize 
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that T c (r) is rather flat over a wide range of r < lO -1 , 
so that any inaccuracy in determining the value of r has 
little effect on T c . 

Our result is in much better agreement with experi- 
ments than T~ , particularly given that hot spot calcu- 
lations tend to somewhat overestimate the value of T c 
compared to full scale calculation without expanding the 
fermionic dispersion near the hot spots^. 

In Fig. [9] we present the eigenfunction of the linearized 
gap equation for various Matsubara frequencies. It shows 
that the infinitesimally small d— wave gap at T c is non- 
monotonic along the FS, with its maximum not at the 
hot spot (the point k y = 0) but at a point away from it. 
This agrees with the result of KC. 



B. The non-linear gap equation 



gular dependence of the gap along the FS, and to verify 
whether it remains non-monotonic below T c . In the pro- 
cess, we also obtain the frequency dependence of the gap 
and compute the ratio 2A max /T c . 

The computational procedure is straighforward. We 
introduce A(k,u n ) and fermionic Z{k,uj n ) via 



T,(k,u} n ) = -iw„(l - Z(k,uj n )) 



A(k,u n ) = Yl 2{k,u) n )/Z(k,w n ). 



(34) 



(35) 



We next derive and solve the non-linear gap equation 
for T < T c . The key goals here are to obtain the an- 



The equations for the normal and anomalous, fermionic 
self-energies are re-expressed via A and Z as 



J 



d 2 k 



UJ m Z(k,u) m ) + e k 



(2tt) 2 ( U 2 m + A{kj Um )2) Z (k, w m ) 2 + e 2 



(36a) 



E 02 (p) = A(p)Z(p) = 
d 2 k 



A(k,oj m )Z(k,uj r , 



(2^) 2 (ojI + A(k, u m ) 2 )Z(k, u m y + e 2 



x(k-p,Lo m - u n ) 



(36b) 



r 



where p = (p,ui n ). We restrict momenta in A to the and we introduced dimensionless k = k y /qo, p = p y /qo, 

and fltot = (Xo/9o) n tot = (4Xo/?o) n Tne lul1 polariza- 
tion operator depends on A(fc,w n ) and is given by 



Fermi surface and introduce 

A(p,w m )| = A(p y ,u m ) 



(37) 



Carrying out the integration w.r.t. as before, and 
using Eq. (|36al) to eliminate Z in the l.h.s. of Eq. ()36bl) . 
we obtain the non-linear gap equation in the form 



A{p,u n ) 



3T 



E 



i 



JST 

4ili, 



Al x(p,k,m n - uj m ) - 

dfc ; =A{k,UJ m ) 

uj 2 n + A 2 {k,uj rn ) 
dfc 2m + 1 Wm ~ 



2n+ 1 



(u4 + A 2 (fc, Wm )) 



where we have 
X(p,k,n n ) 



((p - k) 2 + r 2 (p 2 + k 2 ) 2 + UtotiQ, «n) 



(38) 



(39) 



IItot,sc((9, tin) — 



2T 



E 



where 



dk(E 1 +E 2 + E 3 ) (40) 
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E 1 = -i\l 



E*=i\l 



< + A 2 (fc, Wm ) / (8k 2 - + A 2 (fc, iu m )) 2 + (( Wrn + r>„) 2 + A 2 (fc, iu m + n n )) 



8fc 2 + i w m + J(w m + ft„) 2 + A 2 (fc, w m + fi n ) 



'(w TO + n n ) 2 + A 2 (k,Lo m + Q n ) ) (8fc 2 + (w m + ft„) 2 + A 2 (/c,w m + ft„)) 2 + (cj 2 „ + A 2 (fc,w m )) 

/ 

E 3 = A(fc,w m )A(fc,cj m + fi n ) 



1 



u 2 n + A 2 (k,u m ) ( (8fc 2 - i x /^ 2 „ + A 2 (fc,w m )) 2 + ((w m + a,) 2 + A 2 (fc, Wm + H n )) 



^ (w m + ft n ) 2 + A 2 (fc,u; m + fi„) ^(8A: 2 + iyj (uj m + tt n ) 2 + A 2 {k, u m + ft„)) 2 + (oj 2 n + A 2 (fc, u m )U f 



(41) 



Eqs. (|58|> ) and ([3D"])) have to be solved self-consistently 
for T < T c , to obtain the superconducting gap A(kp, ui n ) 
along the Fermi surface. 



optics we have^i 2A max /T c ~ 5 for Pr 1S5 Ceo.i5Cu04, 
and from Raman measurements^ on Pr2_;cCe x Cu04 and 
Nd 2 - a: Ce a ;Cu04 we have 2A max /T c ~ 3.5. 



C. Numerical solution of the non-linear gap 
equation 

We numerically solve the above set of equations using 
the Gauss-Legendre quadrature. We choose a particular 
value of r ~ 0.08 for which T c ~ 0.10 Q ~ 0.30 u . 

In Fig. [10] we present our results for the solution of 
the coupled non- linear equations (|3"5|) and (|4T)|) . for var- 
ious T < T c . We see that the gap continues to be non- 
monotonic along the FS for all Matsubara frequencies, 
with its magnitude increasing as T decreases, as indeed 
it should. The position of the maximum remains essen- 
tially intact. We consider our results as proof that the 
gap in electron-doped cuprates is indeed non-monotonic 
along the FS. We emphasize that the position of the max- 
imum is not at hot spots, which in our approximation are 
located right where the zone diagonals intersect the FS. 

In Fig. [11] we plot the temperature dependence of the 
gap maximum, at different Matsubara frequencies, mea- 
sured in units of T c . We see that for the lowest Matsubara 
frequency (n = 0) the ratio 

2A max (T^0,n = 0) ^ 4 
T c 

The n = value of A(fcp-,w„) is close to the real fre- 
quency uj = value of the gap; so we predict that the 
ratio of the measured largest 2A along the FS and T c 
should be close to four. This is not far from the BCS re- 
sult for a d-w&ve superconductor—. Experimentally, from 



V. SUMMARY DISCUSSION 

To summarize, we have considered in this paper nor- 
mal state properties, the pairing instability temperature, 
and the structure of the pairing gap in electron-doped 
cuprates. We assumed that the pairing is mediated by 
low-energy, collective spin excitations, and that antifer- 
romagnetic order develops close to the doping where the 
FS touches the magnetic Brillouin zone boundary at four 
points, hot spots, which at this doping lie on the zone di- 
agonals (i.e., antiferromagnetism emerges together with 
the appearance of hot spots). 

Because of the absence of a natural ojjj / Ep parame- 
ter, we extended the low-energy spin-fermion model, us- 
ing Eliashberg theory, to N » 1 fermionic flavors with 
l/N a small parameter. Using this we justified the ne- 
glect of vertex corrections and of dS(k) / de^; and also 
argued for the necessity to solve self-consistently, in the 
normal state, the coupled integral equations for fermionic 
E(cj) and bosonic H(Q,fl). This is necessary because 
the self-energy corrections to the polarization bubble are 
not small in l/N. In this respect, Eliashberg theory in 
electron-doped cuprates is different from the one in hole- 
doped cuprates. In hole-doped systems, the hot spots 
are far from the zone diagonals, Fermi velocities at hot 
spots separated by Q are directed at a finite angle w.r.t. 
each other, and the polarization bubble, II(Q, Q) oc |fi|, 
is the convolution of local Green's functions which do 
not depend on the self energy S(w). In electron-doped 
cuprates, the hot spots at QCP are on the zone diogonals, 
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FIG. 10: The superconducting gap A(k y ,uj n ) (in units of ljo) obtained as a solution of the non-linear gap equation for different 
temperatures T < Tc. We set r ~ 0.08 (T c ~ 0.30 ujo) The index n in the legend correspond to different Matsubara frequencies 
Wn. 



n=0 




0.4 0.5 0.6 0.7 0.8 0.9 1 

T 
T, 



FIG. 11: The momentum peak of A(k v ,u> n ) in units of T c , 
plotted against T/T c . Parameters are the same as in Fig. 1101 



the velocities of fermions separated by Q are antiparallel, 
and the polarization bubble II(Q, Q) oc \/\Q\ is not ex- 
pressed as a convolution of two local Green's functions. 
As a consequence, the self-energy corrections become rel- 
evant, necessitating a self-consistent analysis of H{Q, f2) 
and £(a;). 

In the earlier study by these self-energy correc- 

tions to the bubble were not included, and the T c calcula- 
tion was carried out without taking this self-consistency 
into account. In our work we have computed E and II 
self-consistently. We have argued that, to a reasonable 
accuracy, the self-consistent solution for E and II reduces 
to the renormalization of earlier results by constant fac- 
tors, i.e., II — > all and E — » Ti/y/a, where a ~ 0.6. 

Using these results, we re-derived the formula for T c 
and obtained T c -> T c /a 2 ~ 3T C . Using the same param- 
eters as in the earlier study, we obtained T c ~ 30 — 35K, 
in very good agreement with the experimental data for 
electron-doped cuprates. 

We then derived a coupled set of non-linear equations, 
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for T < T c , for the frequency and momentum depen- 
dent d— wave superconducting gap A(kp,uj n ) along the 
FS and the polarization operator H(Q, These equa- 
tions are highly non-linear as A holds non-linear de- 
pendence on II, and II a holds non-linear dependence 
on A. Our numerical solution of this set of non- linear 
equations shows that A(/cf, w n ) decreases monotonically, 
as expected, with increasing Matsubara frequency w„; 
and more interestingly, that for each w„ the gap is non- 
monotonic in momentum along the FS, with its node 
at the zone diagonal (where a hot spot is also located), 
and the maximum some distance away from the zone di- 
agonal. This non-monotonicity of the superconducting 
gap has been conjectured based on the analysis of Ra- 
man data^ and subsequently detected in ARPES mea- 
surements^ 

For 2A max (T — > 0, n = 0)/T c we obtained a value close 
to four. This is also consistent with the data: from op- 



J 



tical measurements 2A max /T c ~ 
Raman measurements 2A max /T c 



5 (Rcf. 31 ), and from 
-3.5 (RefiS). 
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Appendix A: Calculation of II V 

In this appendix we present the details of the calcula- 
tion of the polarization operator with a vertex correction, 
Eq. d5IJ). We have 



n v (Q,n) = -2Ng 4 



00 dwd 2 fcdu/d 2 fc' 



x(k — k' ' ,uj — co') 



(2tt) 6 (iw - e fe )(i(w + O) - e k+Q )(iu' - e fe <)(i(a/ + Q)- e k , +Q ) 



(Al) 



Using Eqs. ([3]) and (Q} for e k and e k +Q and carrying out 
a change of variables, we re- write (|A1[) as 

I 



n Y (Q,n) 



2Ng* 



dcoduj'dky dk' y de^de^ 



x(k' — k, lj' — oj) 



(e fc - iw)(e fc - 2/3 2 fc 2 + i(w + Sl))(e k . - iw')(e k . - 2j3 2 k' 2 + i(w' + (I)) 



(A2) 



r 



Once again, neglecting the dependence of x on e k , we set, 
X(k' -k,w' -u)) = x(k' - k,u' - w)| £fc=efc , =0 (A3) 

Using equation © for %, and eliminating the dependence 
on k x , k' x using 



£fc = v F k x +/3 2 k 2 = 



v 

e k > = -vfK + /3 2 k' 2 = 



(A4) 



we obtain, 



X{k' y ,ky,u)' -u) = x{k' -k,u}'- w)| £fc=£fc , =0 

= Xo 

(ky - k' y Y + + fc; 2 ) 2 + 4xo n(Q, u > - u) 

Xo 



(A5) 



q 2 ( (k - k 1 ) 2 + r 2 (fc 2 + fc' 2 ) 2 + jV^/i^-M 



In the last line we have replaced k y and k' y by dimension- 
less k = ky/qo and fc' = k' y /q n , and used Eq. ([IT)]) for II 
with Fn — 1 because we are at T — > 0. Also the ujq that 
appears is as given by equation 0, i.e. without the large 
iV redefinition. We have verified that for the present cal- 
culation the r dependent term in x is irrelevant and can 
be dropped. The integrations w.r.t. e k , e k > are factorized 
and easy to carry out. The result is 
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n v (Q,J2) 



v 2 F (2nf 



2th 



dkydk' y 



duj' 



2-ki 



o i(2w + O) - 2/3 2 /c 2 



2?ri 



da; 



o i(2w' + n) - 2/3 2 fc' 2 



2?ri 



-n 



doj' 



2$ 2 k 2 - i(2a/ + Q) 



2fi 2 k 2 y - i(2w + 0) 
X(fcy - k y ,u}'-u) 



(A6) 



Re- arranging the integrals over dw and da/, we get a more compact result 



n v (Q,n) 



27V ff 4 



dfcj,dfcy( / dwdu/ 



dw duj'){D(uj,k y )D{uj , ,k' y ) + c.c.)x{k' y -ky,w'-w) (A7) 



where 



and 



D{uj,k y )D(cu',k y ) +c.c. = 2 



2/3 2 fc 2 - i(2w + fi) 
2p 2 kpP 2 k' y 2 - (2cj + n)(2o/ + n) 

((2/3 2 fc 2 ) 2 + ( 2w + r ,)2 ) (( 2; g2 fc ,2 ) 2 + (2a/ + n) 2) 



(A8) 



(A9) 



Substituting this into (|A7|) and using 2(3 2 k 2 = 8wo£; 2 , where fc is the dimensionless variable introduced above, we 
obtain 



n v (Q,n) 



d/cdfc' 



F n (fc,fc',w,w / ) 



(A10) 



(fc - A;') 2 + JVVIw- w'l/wo 



where 



F a (fc,fc',w, w ') 



( fc4 + (^) 2 )( fc ' 4 + (^) 2 ) 



(All) 



Extending this result to large N as described in the main text (Eq. (|18p). i.e., redefining u!q — > ujq/N 2 ; and 
introducing new variables 



and 



we obtain 



x = uj/oj q , y = uj'/uj 



k = Nk, k' = Nk' 



(A12) 



(A13) 



n v (Q,o) 



2tt/3v f 64tt 2 
dfcdfc' 



( / dxdy dx dy) 

Jo J — oc 

Pfc' 2 - AS 



(A14) 



(¥ + A 2 )(k> + B 2 )((k - k') 2 + 2N 2 y/\A-B\) 



where 



.4 
B 



2x + W 
8 

2y + il' 



(A15) 
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and 



Rearranging the limits of integration in the second set of x, y integrations, we obtain, 



(A16) 



n v (Q,n) 



N 2 g 2 
2tt/3vf Mir 2 J 



dxdy / dfcdfc' 



fc 2 fc' - AB 



(fc 4 + A 2 )(k' + B 2 )((k - fc') 2 + 2N 2 V / \A-B\) 



k 2 k' + AB 



(A17) 



(fc 4 + A 2 )(k' + B 2 )((k - fc') 2 + 2N 2 y/\A + B\) 
Performing the integration over fc and fc' we obtain, after some algebra, 



n v (Q,fi) 



2tt/3w f 64 



J dxdy ( 



I VI- VB| 



\ VAB(AN 4 (y/A + y/B) + \VA- y/B\ 3 ) 

A- B + 2N 2 y/\A- B\ 



2NVBJ V / \A -B\ \ [A-B- 2NVBJ V / \A - B\ + 4N 2 y/\A - B \( N J V\ A ~ B \ + ^B) 2 



B - A + 2N 2 y/\A- B\ 



2Ny/AJ^\A -B\ \ LB - A- 2Ny/J J ^\A - B\ + 4N 2 y/\A - B\(y/A + NJ y/\A - B 



2(VAB + N 2 ^/A + B) 



IAB{{A + B) 2 + 4N 4 (A + B) + 8N 2 VAB^/A + B) 

A + B- 2N 2 y/A + B 



2NVBy/VA + B[(A + B + 2N v r By / ^A + B) + 4N 2 y/A + B(VB + N\/ \J A + B) 2 



A + B- 2N 2 VA + B 



2Ny/AVVA + B I (A + B + 2N^A^J^A + B) + 4N 2 y/A + B{\[~A + N y/VA + B) 2 



(A18) 



We note that our result is symmetric w.r.t. to A <-> B This is the result we quoted in the main text, Eq. (|22l) . 
as it should be, since the original integrals are symmetric 
w.r.t. uj -H- uj' . Evaluating the remaining integrals, sub- 
tracting the Q— independent term coming from large x, y 
(i.e., from high internal frequencies), and taking large N 
limit, we find the leading term to be 
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